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ITRODUCTION 

Delta  modulation  is  a  hybrid  system  composed  of  pull 
duration  and  pulse  amplitude  modulation  in  conjunction  with 
feedback  which  encodes  a  signal  into  a  sequence  of  quantized 
pulses. 

This  report  reviews  previous  work  on  the  subject  which 
ranges  from  descriptive  presentations,  to  probabilistic  analyses 
from  the  information  standpoint,  to  newer  concepts  of  the 
character  of  delta  modulation.   It  then  presents  an  ideal  mathe- 
matical model  for  a  delta  modulation  system  and  a  resulting 
deterministic  method  of  analysis.   The  results  of  this  analysis 
are  given  and  arc  followed  by  a  discussion  of  improvement 
techniques  on  the  basic  system.   A  comparsion  of  the  deterministic 
and  probabilistic  approaches  is  made,  and  it  is  shown  that  the 
two  approaches  are  complementary  in  scope.   Finally  a  generali- 
zation of  delta  modulation  to  a  system  having  five-level  output 
quantization  is  made. 

PREVIOUS  WORK 

F.  de  Jager's  Paper 

F.  de  Jager  (4)  gave  a  qualitative  description  of  the  delta 
modulation  system.   No  extensive  mathematical  analysis  was 
attempted. 

The  physical  system  (sec  Fig.  1)  consists  of  a  comparator, 
D,  which  controls  a  pulse  modulator,  PMt  fed  by  a  pulse  generator, 
PG.   The  signal  from  the  modulator  output  passes  through  an 


integrating  network,  F^ ,  and  back  to  the  negative  side  of  the 
comparator.   The  delta  modulated  signal  is  taken  from  the 
pulse  modulator  output,  and  after  passage  through  a  channel,  the 
original  signal  is  reconstructed  by  means  of  an  integrating  net- 
work, F2,  similar  to  F^. 

de  Jager  found  that  using  a  single  integrating  network  for 
F^  gave  satisfactory  results.   In  an  effort  to  improve  the 
signal-to-noise  ratio  he  tried  using  a  second  integrator  in 
cascade  with  the  first.   This  resulted  in  noise  reduction  but 
gave  rise  to  stability  problems  and  difficulties  in  transmitting 
the  higher  frequencies.   He  compromised  by  using  a  single  inte- 
grator followed  by  a  low  pass  filter  with  cutoff  at  3800  cps. 
Actually,  the  single  integrator  was  a  shallow-skirted  low  pass 
filter  with  cutoff  at  200  cps. 

The  difference  between  the  original  and  the  reconstructed 

signals  was  defined  as  quantizing  noise.   Using  this  definition 

signal-to-noise  ratio  formulas  were  calculated  for  both  single 

and  double  integration.   For  single  integration, 

3       1 
S/M  =  Cx   fj7  /  f  f02  (1) 

For  double  integration, 

S/N  =  C2  fiZ  /  f  f0  (2) 

where   f   is  the  frequency  of  the  input ,  f ^   is  the  pulse 
frequency,  and   fQ   is  the  cutol ;      ;uency  of  the  low  pas6 
filter.   The  C^'s  must  be  determined  by  numerical  methods.   A 
comparison  of  these  theoretical  rv  Milts  with  experimental  values 


showed  that  the  improvement  in  signal  to  noise  ratio  as  pulse 
frequency  increased  was  in  fair  agreement  with  the  5/2  and  3/2 
:  oi  the  pulse  frequency  for  double  and  single  i;      tion 
respectively. 

The  necessary  bandwidth  for  delta  modulation  was  stated  to 
be  about  50%  greater  than  for  pulse  code  modulation.    The  quality 
of  DM  was  said  to  be  roughly  the  same  as  that  of  an  eight-unit 
code  PCM. 

L«  h.  Zetterberg's  I aper 

L.  H.  Zetterberg  (9)  compared  delta  modulation  and  pulse  code 
modulation  from  the  point  of  view  of  principle  as  well  as 
performance,   lie  concluded  that  they  are  similar  because  both 
represent  a  signal  by  means  of  a  sequence  of  pulses  and  further- 
more because  the  systems  may  be  characterized  by  the  same 
parameters.   This  comparison  may  not  have  been  completely  valid, 
but  at  least  PCM  served  as  a  standard,  and  Zetterberg  was  able 
to  obtain  a  measure  of  the  performance  of  DM. 

In  accordance  with  Shannon's  general  ro.tnod  of  analyzing  a 
communication  system,  delta  modulation  was  viewed  as  a  discrete, 
noiseless  channel  and  channel  capacity  was  calculated  using  this 
assumption.   The  source  of  information  was  considered  to  be  a 
set  of  functions  generated  by  finite  Markov  chains.   A  Markov 
chain  is  a  stochastic  process  in  which  a  future  state  is  deter- 
mined only  l^y  the  previous  state  and  the  transition  probabilities 
remaining  fixed.   For  speech  signals,  or  any  signal  likely  to  be 
transmitted,  there  will  be  a  much  stronger  dependence  between 


the  various  states.   But,  according  to  Shannon,  a  stronger 
dependence  between  the  values  of  the  functions  gives  a  smaller 
entropy  for  the  process  and  thereby  a  lower  rate  of  information. 
Therefore,  use  of  the  Markov  chain  will  give  only  an  upper 
bound  on  the  channel  capacity ,  but  even  this  much  will  be  very 

useful. 

Zetterberg's  results  showed  that 

G  =  i  log  2  (3) 


for  PCM,  and 


C  =  |  log  (2  cos  fr/n+1)  (4) 


for  DM,  where  G  =  channel  capacity,  T  =  pulse  interval,  and 
n  =  number  of  quantizing  levels,  i.e.,  the  quantized  output 
magnitudes. 

It  can  be  observed  that  as  n  increases  to  10  or  greater,  the 
channel  capacity  of  a  DM  system  approaches  that  of  a  PCM  system 
so  closely  that  the  difference  is  almost  negligible.   Since  in 
a  practical  application  the  number  of  levels  is  likely  to  be 
much  greater  than  ten,  deltA  modulation  can  be  seen  to  compare 
favorably  with  regard  to  channel  capacity. 

The  paper  continued  with  spectral  calculations  from  which 
few  conclusions  can  be  drawn.   I'hen  analyses  of  the  error  signal 
were  undertaken.   Two  major  causes  of  error  are  the  sampling 
process  and  the  quantization  of  the  output.   Inherent ly  sampling 

•ulta  in  error  which  can  bt  reduced  mainly  by  increasing  the 
sampling  rate.  Zetterbcrg  divided  the  eit'ects  of  quint izat ion 
into  those  caused  by  overloading  and  the       »>d  by  granulation. 


For  DM,  overloading      rs  when  the  slope  of  the  input  signal 
is      *r  than  tbi  height  of;  one  DM  increment  per  unit 
sampling  time.   For  PCM,  overloading  appears  when  the  amplitude 
ot  the  signal  becomes  juv.iLcr  than  ttv     imum  of  the  quantizing 
code.   However,  Zctterbcrg  failed  to  point  out  directly  that 
essentially  the  same  effect  occurs  in  DM  when  the  input  amplitude 
becomes  greater  than  the  saturation  level.   The  other  effect, 
nulation,  is  an  interpolation  error. 

Values  of  noise  power  were  calculated  for  both  overloading 
and  granulation  effects.   Signal-to-noisc  ratios  were  then 
computed  for  the  various  cases  .->.nd  the  results  for  DM  were 
compared  with  those  for  PCM.   It  was  found  that  in  certain  cases 
the  noise  power  due  to  granulation  errors  is  very  much  less  with 
delta  modulation  than  with  PCM.   Furthermore,  for  higher 
frequencies,  overloading  in  DM  will  be  less  troublesome  than  in 
::,  and  it  occurs  much  more  quickly  in  PCM. 

Considering  bandwidth  requirements,  it  was  found  that  the 
channel  bandwidth  for  DM  must  be  approximately  double  that  for 
rCM. 

B.  L.  Barber's  Iteport 

Barber's  master's  report  (1)  gave  different  insipht  into 
the  n.iture  of  delta  modulation.   This  report  established  the 
basis  for  a  deterministic  analysis.   Previous  writers  have 
compared  DM  to  PCM  using  the  latter  for  a  yardstick.   Although 
they  did  not  represent  DM  as  actually  being  a  form  of  PCM,  they 
maintained  that  the  systems  are  quite  "similar".   However,  Barber 


has  shown  DM  in  reality  to  be  derived  from  pulse  duration 
modulation  and  pulse  amplitude  modulation. 

The  pulse  duration  modulator  is  constructed  as  shown  in 
Fig.  2.  The  input  signal  is  added  to  a  triangular  wave  and  the 
sum  is  fed  into  a  signum  function,  whose  output  is  the  PD>: 
signal.   Demodulation  consists  of  averaging  the  PDM  signal  over 
the  period  of  the  triangular  wave.   It  can  be  shown  that  this 
average  is  proportional  to  the  modulator  input  level  up  to  a 
saturation  level.   As  will  be  shown  later,  such  a  demodulator 
can  be  realized  by  using  a  holding  filter  with  holding  time  equal 
to  the  period  of  the  triangular  wf.ve. 

Pulse  amplitude  modulation  is  accomplished  by  multiplicative 
sampling  of  the  input  by  means  of  a  periodic  delta  function. 
A  PAM  demodulator  is  ideally  a  steeply-skirted  low  pass  filter. 
If  the  period  of  the  sampling  function  is  T,  the  cut-off  frequency 
of  the  low  pass  filter  will  be  l/r. 

Barber  proposed  to  let  the  PDM  and  its  demodulator  be 
separated  by  a  PAM  and  its  demodulator.   (See  Fig.  3.)   He  called 
this  combination  of  PDM  nnd  PAM  "pulse  duration-amplitude 
modulation".   The  new  PDAM  modulator  consists  of  the  cascaded  PDM 
and  PAM  modulators,  and  the  new  demodulator  is  a  low  pass  filter 
in  cascade  with  a  holdin   filter. 

If  the  triangular  w  ive  is  removed  from  I        modulator 
the  weakened  modulator  can  be  soneuh.it  improved  by  negnti- 
feedback  so  that  the  output  of  the  holding  filti     places  t 
triangular  wave.   This  arrangement  is  shown  in  Pi|  .  •* .   It  can  be 
seen  that  this  hybrid  sysU-m  is  fftlaai    -  tu  mo     urn. 
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Other  Papers 

E.  E.  David  (3)  developed  interesting  results  from  studies 
of  delta  modulated  encodings  for  transmission  of  pictures.   He 
also  called  this  "differential  quantization".   The  quantizer 
used  for  the  picture  encoding  was  eight-leveled  with  the  levels 
tapered  and  equivalent  to  3  bits  per  picture  point.   The 
performance  of  this  was  compared  to  that  of  a  3-bit  pulse-code 
modulation.   Delta  modulation  produced  a  superior  reproduction 
of  the  picture.   PCM  reproduction  had  objectionable,  ragged 
contours  in  the  picture.   Here  is  a  case  where  the  quality  of 
delta  modulation  outstrips  that  of  its  standard,  pulse  code 
modulation. 

A  paper  by  Lender  and  Kozuch  (6)  provided  experimental 
data  on  the  operation  of  delta  modulation.   A  discussion  of  the 
hardware  used  in  the  construction  of  an  experimental  model  was 
given,  and  an  evaluation  of  the  performance  was  made.   This 
evaluation  was  made  by  means  of  signnl-to-noise  ratios  for 
sinusoidal  inputs  and  by  means  of  subjective  tests  with  listeners 
for  voice  inputs.   These  results  were  compared  to  similar  ones 
pulse  code  modulation  over  n  ran"*-  of  sampling  rates.   The 
cor;  DWad    at  delta  modulation  is  bettor  th.     H  for 

impllng  rates  of  40,n  oiul  or  less  ami  that  above 

this  rat>    M  is  bettor.   Ho     :,  Li  bandwidth  is  no  problem, 
DM  can  be  niadc  to  match  7-bit  i  CM  by  Uiii  width  40fc 

greater  than  the  bandwidth  require!  tor  tlu* 
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Fig.    Ij..       Delta   modulation   as   derived   from  PDAM. 
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DETERMINISTIC  ANALYSIS 


Z-transform  Description  of  Modulator- 
Demodulator  Model 


The  description  of  the  delta  modulation  system  given  by 
Barber  provides  a  basis  for  devising  a  mathematical  model  in 
terms  of  Z-transforms  by  means  of  which  one  can  make  a  deter- 
ministic analysis.   The  same  degenerate  pulse  duration  modulator 
may  be  used  for  a  model.   The  pulse  amplitude  modulator  is  the 
sampling  operator,  Z. 

The  PDM  demodulator  is  easy  to  represent.   It  has  been 
mentioned  that  the  requirement  for  this  demodulator  is  an 
integrating  network  which  averages  over  the  period  of  the  tri- 
angular wave.   Let  this  period  be  mT,  where  m  is  an  integer  and 
T  is  the  period  of  the  sampling  operator.   If  g(t)  is  the  function 
to  be  averaged  over  mT ,  then  the  average  at  time  t  is 

— T7T-       g(t)  dC  .   Assuming  that  g(t)  is  piecewise  continuous, 

^t-raT 

taking  the  Laplnce  transform  of  the  averaged  function  gives 

-fe  (L  [tft>]}  (I-.-*-)      -     izg-S     L[,<t>]  (5) 

,    ~-mTs  ,    ^.m 

Therefore   the    required  operator   is      -=^ =     — - —    ,    U'hore 

mTs       mTs   " 

z  =  e~  a),  which  is  the  transfer  function  of  a  holding  filter. 

The  PAM  demodulator,  a  steeply-skirted  low  pass  filter,  is 
not  satisfactory  for  stable  o;     ion  of  a  feedback  device,  and 
a  filter  with  shallow  skirts  must  be  used.   This  simply  degrades 
fchi  Original  com  i      ion  a        rther.   Another  holding  filt 
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Cutoff  frequency  l/T  can  be  substituted.   Its  transfer 
function  is  (l-e"Ts)/s  =  (l-z)/s.   To  show  that  this  is  a 

illow-skirtod  low  pass  filter,  replace  s  by  jW  in  the  transfer 
function  and  find  the  spectrum. 

4  ■    »    l-%"         .    TW  m      sin  Q 

T(JW)  .  — ^—  -  -  .in  -  .  t  —£—  (o) 

(where  0  =  ~t)»   This  function  has  a  shallow,  low  pass 
characteristic  with  cutoff  occurring  at  9  -  TT  or  f   =  l/T. 
The  phase  angle,  arg  T(jW)   =  W/2 ,  is  linear,  assuring  a 
constant  delay  for  all  frequencies.   The  PDM  holding  filter  has 
similar  response  characteristics  with  a  cutoff  frequency 
fc  =  1/mT. 

The  Z-transform  model  of  the  delta  modulation  system  now 
appears  as  shown  in  Fig.  5,  and  is  in  suitable  form  for  analysis. 

Recursion  Formula  for  Model's  Output 

The  objective  of  a  deterministic  analysis  is  to  determine 
the  output  of  the  demodulator  as  a  function  of  time,  y(t),  given 
the  modulator  input,  x(t).   A  solution  yielding  the  complete 
function,  y(t),  will  be  very  difficult  to  obtain,  but  by  means 
of  z-transform  methods  it  is  possible  to  find  the  value  of  y(t) 
at  each  of  the  sampling  times,  t  =  nT,  where  n  =  1,  2,  .  .  .  . 
For  notational  purposes  let  the  Laplace  transforms  of  x  and  y 
be  x  and  y   respectively.   The  derivation  proceeds  as  follows: 
the  input  to  the  signum  function  is   (x*  -  y) .   Then 


1-z 

1-z 
s 

mTs 
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y  =[  kZ   sgn  (x  -  y)J        (7) 

See  the  definition  of  sgn  £"  on  page  13. 

Now  take  the  Z-transf orra  of  both  sides  of  the  above  relation, 
observing  that  Z  may  be  taken  inside  the  signum  function  since 
the  order  of  sampling  and  quantization  is  immaterial. 


Zy  =  mT  Z fe  [(l-zm)(l-z)  sgn  (Zx  -  Zy)] 

By  use  of  the  relation, 

Z  |f(s)  Z  [g(s)]J  =  Z  [PCs)]  z[g(s)] 
it  follows  that 

(l-zm)(l-z)  sgn  (Zx  -  Zy) 


k 

^  "  mT 


z  (±z) 


; 


:iow  z(s^)=H7f^  • 


(8) 


(9) 


(10) 


Then 


or 


Zy  =  m  (X^-Jz  sgn  (Zx  -  Zy) 


(11) 


Zy  -  zZy  =  H  [z  sgn  (Zx  -  Zy)  -  z™*1  sgn  (Zx  -  Zy)] 

(12) 

But  since  z  is  a  delay  operator  it  may  be  taken  inside  the 

signum  function  giving 

Zy  -  zZy  =  m  [svm  (zZx  -  zZy")  -  sgn  tz"1*1  Zx  -  zm*1    Zy)J 

(13) 

By  writing  this   equation  in   terms   of   the    sampled   and  delayed 

values    of  x   and  y,    one   obtains    the   recursion  ion, 

Vn    "   Vl   =   BP*11    (xn-l    ■   f^)    ■    SZ*  -1    ■    :■.._..- 1>] 

(14) 
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Discussion  of  Results  for  Step  Inputs 

Co'   .•  r  i'ut.  put  s .   The  previous  recursion  formula  is 
particularly  well  suited  for  computation  on  a  digital  computer. 
It  is  Mfc<  •    ry  only  to  specify  the  values  of  xu  and  the  value 
of  yQt  and  the  computer  will  do  the  rest. 

This  particular  set  of  calculations  w     rri<«d  out  on  an 
IDK  650  digital  computer.   rrovision  was  made  in  the  programming 
to  compute  data  for  various  types  of  input  functions  and  for 
different  values  of  system  parameters.   The  answers  were  supplied 
by  the  machine  at  the  rate  of  approximately  500  per  minute. 

The  following  system  parameters  were  used: 
=  2,   m  =  256,   T  =  1 
The  signum  function  was  such  that 

f  i ,  £ > o] 

sgn?=/  0,  5  =  0  (15) 

This  does  not.  agree  completely  with  the  models  of  de  Jager  and 
Zetterberg,  which  were  binary  in  character.   They  were  such  that 
the  output  for  zero  input  was  a  series  of  alternating  positive 
and  negative  pulses.   Such  a  characteristic  could  be  obtained 
here  by  assigning  either  +1  or  -1  to  sgn^  when  ^  =  0.   The 
performance  of  a  three-level  device  can  be  expected  to  be 
somewhat  better  than  for  the  two-lev<;l  case. 

Results  of  the  calculations  showed  the  overall  input- 
output  characteristic  to  be  almost  that  of  a  servo  for  values 
of  x  between  k/m  and  k.   Below  x  =  k/m  there  was  no  response. 
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Above  x  =  k  saturation  occurred.   In  the  useful  range  the  out- 
put, after  reaching  the  input  level,  contained  fluctuations. 
For  the  particular  case  calculated,  y  built  up  linearly  to  x  and 
remained  equal  to  x  until  sample  number  257.   Thereafter  it 

dropped  to  (x  -  — )  and  remained  at  that  value  for  a  number  of 

ra 

intervals  depending  on  the  value  of  x.   It  then  returned  to  x 
and  remained  until  sample  number  514  where  it  again  dropped. 
This  pattern  appeared  to  repeat  indefinitely.   See  Table  1  for 
rise  times  and  durations  of  the  drops  (or  holes)  for  different 
inputs,  and  see  Fig.  6  for  a  typical  curve.   Since  the  period 
of  the  degenerate  triangular  wave  was  raT ,  it  can  be  seen  that  the 
period  of  the  sets  of  holes  is  about  the  same  as  that  of  the 
triangular  wave.   The  reason  for  the  fluctuation  could  be 
explained  as  the  need  tor  the  feedback  system  to  generate  its 
own  "triangular  wave". 

Errors  and  Signal-to-Noise  Ratios.   For  inputs,  x,  restricted 
to  powers  of  2   an  empirical  relation  between  the  duration  of 
the  drop  and  the  input  magnitude  can  be  made.   From  the  tabu- 
lation it  can  be  seen  that  the  number  of  drop  intervals  is  an 
integral  power  of  2.   Then  if  x  =  l/2n  for  k/m  <:   x  <  k,  then  is  , 
the  number  of  drop  intervals,  is  equal  to  2  ~n.   For  example, 
for  x  =  1/4 ,  n  =  2  and  v  -  2   "'  =  32,  which  is  in  agreement  with 
the  tabulated  values.   Now  the  magnitude  of  y  is  ~  -  -  =  - — "  T2U 
for  2  ~n  intervals  and  Tfi  for  the  remaining  257  -  27"n  intervals 
in  each  period.   Hence  the  aver.!.*-  value  of  y  can  be  compute 
as  follows: 
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- 

"  r  l  ^  -  2',<)  rial 

7  (  Jn'  ~  Tn   x  (I6) 

>n  this     1  n  l»e       that  thi        e  value  of  the  OU 

lid  be  r  *  <V  equal  exactly  to  the  value  of  the  input  by  means 

of  a  suitable  gain,  attaining  a  servo  response.   (This  has  been 

shown  only  for  values  of  x  equal  to  2n.) 

Subtracting  the  average  signal  from  the  total  eigtial 

yields  a  noise  signal  which  will  enable  one  to  compute  a  signal- 

to-noisc  ratio.   Let  the  amplitude  of  the  positive        D  of 

the  noise  signal  l>e  called  u.   Then  u  ■  x  -  (y)  =  -s^  -  r^r=  •  5TT  a 

— =—  .  -=—  .  The  amplitude  of  the  negative  portion  will  be 

25  7  2n 

1  o7-n 

T28~  -  u  =  Y23  ^  "  257  ^*  Tnen  tne  ™s   vaiue  of  tne  noise 

signal,  N,  is  given  by  the  following  equation: 


V^fljw  -  ^T-^^-hY  (-7  -^1    (w 


Simplifying  yields 


5 


■^f\j?.:.7-2"-7  -l'  (18) 


The  signal  to  noise  ratio  is  then 

_   256 

s/N  =v277~2n-7  .r  (19) 

>ressed  in  decibels  this  is 

S/N  (db)  =  48.16  -  10  log10  (257- 2n~7  -  1)       (20) 

One  must  remember  that  the  preceding  are  empirical  relationships 
used  to  cxi/cdite  noise  calculations.   See  Fig.  7  for  a  curve 

of  calculated  signal  to  noise  ratios  plotted  versus  step  input 
magnitude  . 
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Table  1.   Rise  tine  and  number  of  hole  Intervals 
for  step  Inputs,  m  ■  256. 


Inpui 


2 

1 

1/8 

1/16 

1/32 

l/l2$ 
1/256 


Rise  Time 


Humber  of  holes 
per  257  samples 


25#T 

256 

123T 

128 

61|T 

6k 

32T 

1? 

16t 

16 

8t 

8 

l|S 

I* 

2T 

2 

T 

1 

Table  2.   Ri3e  time  and  number  of  hole  intervals 
for  3tep  inputs,  m  n  16. 


Input 

:    Ri 

se  Time 

::umbor  of  holes 

;     per 

17   samples 

2 

16T 

16 

1 

8t 

8 

*# 

l|T 

h 

\% 

2T 
T 

2 

1 

1/16 

- 

- 

1/16 


■\ r 


. — ^ 


t/T 


256  '  /       512 


PiC.  6.   Output,  y(t),  for  atop  input,  x(t)  -  l/l6, 
n  ■  256. 
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short  set  of  calculations  for  step  inputs  with  m  =  16 
and  k  =  2  was  made.   The  result  :l.itcd  in  Table  2, 

showing  rise  times  and  durations  of  drops.   It  can  be  seen  that 

alts  were  very  similar  to  those  for  m  =  256,  with 

corresponding  outputs  having  similar  properties.   The  period 

of  output  oscillation  was  17  s    .  . .   A  more  general  set  of 

empirical  formulas  can  now  be  determined.   The  period  of  output 

oscillation  is  equal  to  (m  ♦  1).   The  duration  of  the  drops  is 

equal  to  xm/2 ,  where  x  and  m  are  both  integral  powers  of  2. 

Formulas  for  the  average  output  magnitude,  the  rms  noise  value, 

and  the  signal-to-noise  ratio  could  be  derived  from  these  more 

general  relations. 

Signum  Function  with  Dead  Space.   A  few  calculations  using 

a  dead  space  of  magnitude  5  ^n  tae  signum  function  were  made. 

This  new  function  was  defined  as 

f   1,    |>6l 
sgn     ?  =    \     0,     |5|  <;  6  (21) 

6       I-1.  T<-6j 

In  the  delta  modulator,  the  old  sgn  £  was  replaced  by  sgn  ^ 
with  no  other  system  changes. 

For  &<k/m  no  change  was  observed  in  response  to  step 
inputs.   For  6=  k/m  the  magnitude  of  y  was  reduced  by  an  amount 
equal  to  k/m.   Furthermore  the  number  of  hole  intervals  per 
(m+1)  samples  was  reduced  by  1 .   As  Tight  be  expected,  y  =  0  for 
X  =  k/m.   These  results  indicated  that  there  is  no  advantage  to 
a  dead  space  alone.   A  modification  of  the  dead  space  to  produce 
a  five-level  signum  function  and  the  necessary  changes  in  the 
system  for  such  a  modification  will  be  explored  later  in  the  paper, 
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Delta   Modulate  |  s    for     it  i  |     In]    .'.  g .      To   r.akc    the 

on  more  nearly  complete,  a  calculation  of  the  modulator 
output  function  was  made.   This  is  known  to  be  a  sequence  of 
positive  and  negative  pulses.   With  the  signum  function  defined 
to  be  equal  to  0  when  its  argument  is  zero,  or  with  a  dead  space, 
the  modulator  output  becomes  a  sequence  of  positive,  negative, 
and  null  pulses. 

To  obtain  a  relation  between  the  modulator  output  and  the 
input  consider  the  modulator  shown  in  Fig.  5. 

Let  k  be  placed  in  the  feedback  return  for  calculation 
convenience.   The  equation  of  this  system  is 


where  ^  is  the  modulator  output  function.   But  this  is 
equivalent  to 

sgn   Zx    &-  (l-Zm)(l-z)  Z(^  $)]       (23) 


i  c 


mi 


(24) 


Now  zf— 7^  lis  difficult  to  express  as  a  function  of  the 
z-transforras  of  — j  and  c  .   But  Halijak  (5)  has  shown  by  his 
integrator  substitution  program  that  such  an  operation  may  be 
approximated  by 

z  [^?  <s)]  ;  T  z&)[zr  (s)  "  ^f  (t) 

and  hence 

Therefore  the  substitution  can  be  made 

^=   Sgn     .Z*   "   mT~   (1-zn)(1-2)   flTzT?      (ZS-^o>] 


t  =  0 


(25) 


(26) 
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Now  by  taking  the  z-transform  of  both  sides  and  observing  that 
£0  =  0,  the  following  relation  results 

Zf  =  sgn[Zx  -  S  (l-z^z  Z^J  (27) 

Z^  =  sgn[zx  -  T\^l   zr)z'?]  (28) 

Writing  this  equation  in  terms  of  the  sampled  values  of  x  and  *- 
produces  the  recursion  formula. 

5."  «■»[*- *£$■*]  (29> 

This  recurrence  relationship  was  programmed  for  the 
IBM  650,  and  modulator  outputs  were  computed  for  several  values 
of  x,  with  and  without  dead  spaces  in  the  signum  function. 
Calculation  shows  that  a  positive  pulse  is  emitted  for  each 
interval  until  y  equals  x,  on  the  initial  build-up.   Then  no 
pulses  are  emitted  until  sample  238  when  a  sequence  of  consecutive 
pulses  is  emitted.   The  number  of  pulses  in  this  sequence  turns 
out  to  be  exactly  equal  to  the  number  of  hole  intervals  in  the 
output  y.   This  sequence  of  positive  pulses  is  followed  by 
another  sequence  of  0  pulses.   The  next  pulse  occurs  at  sample 
515  with  the  same  pattern  as  before  repeated.   The  entire  pattern 

repeated  indefinitely  with  ■    riod  of  257  samples.   Con.- 
the  samples  of  yn  which  are  not  equal  to  x.   If  those  I     Laycd 
by  an  amount  equal  to  I ,  they  will  occur  coincidently  with  the 
positive  puis*!  from  the  modulator.   This  seems  rea*.     i.e  when 
one  studies  the  system  and  observes  that  the  future  output  of  the 
signum  function  depends  on  the  presv     .  lue  of  the  error. 


21 


view   :    ..ults  for  Step  Irruts.   A  brief  review  of  the 
deterministic  analys     ows  that  a  mathematical  model  of  the 
system  can  be  construe  rod f  from  which  a  formal  relationship 
between  output  and  input  can  be  derived.   By  programming  a 
. Lai  computer  to  solve  for  the  outputs  corresponding  to 
various  step  inputs,  it  is  possible  to  make  an  evaluation  of 
the  performance  of  the  delta  modulator.   Between  certain  input 
amplitude  limits  the  system  functions  as  a  servo  with  the 
output  equaling  the  input  after  a  finite  delay.   A  measurement 
of  tho  noise  ] ower  generated  in  the  system  is  also  available 
allowing  a  computation  of  signal-to-noisc  ratios  from  the 
various  inputs.   Finally,  the  form  of  the  modulated  signal  is 
computed  by  simple  changes  in  the  formula. 

Performance  and  System  Modifications  for  Ramp  Inputs 

(1-A)+Az  Feedback  Transfer  Function.   Since  the  delta 
modulation  system  can  be  made  to  act  like  a  servo  for  properly 
chosen  step  inputs,  an  obvious  extension  of  the  investigation 
is  to  determine  the  response  for  appropriate  ramp  inputs.   One 
can  be  sure  that  the  slope  of  the  ramp  can  never  exceed  2k/mT 
since  |yn  -  yn_j|<2k/mT.   Also  the  magnitude  of  the  output 
cannot  exceed  k  since  this  is  the  saturation  level  of  delta 
modulation  output.   To  show  this,  consider  equation  (11),  and 
let 

z  sgn  (Zx  -  Zy)  =  sgn  Z  P 
where  ?■  (xn_L  -  y^). 
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Then 

k  i  -7ra      —   k  n~^   >- 

zy-ibrs^z^-^^0z  B^z?  (30) 

In  terms  of  sampled  values  this  is 


k  m"1 


Vn   =  S  ZT  sg*  ?n.; 


(31) 


r=o 

But   -1  ^  sgn^  <  1 

Hence   -k  <  yn  <  k  and  since  y  lies  between  these  bounds  it 
cannot  follow  x  beyond  them. 

It  was  suspected  that  the  output  would  tend  to  lag  the 
input,  and  compensation  for  this  was  attempted  before  any 
calculations  were  made.   In  order  to  bring  the  magnitude  of  the 
output  closer  to  the  inTut ,  the  difference  signal,  £  ,  seen  at 
the  signum  function  could  be  increased  or  be  made  positive  more 
frequently.   Since  ramps  are  the  inputs,  the  introduction  of  a 
delay  into  the  loop  return  would  result  in  earlier  values  of  y 
being  subtracted  from  x,  and  thus  increasing  the  average  of  £. 
Such  a  delay  can  be  introduced  by  placinn  the  delay  operator,  z, 
in  the  loop  return.  P"n  would  then  become  (xn_^  -  vn-2)»   ^en 
better  might  be  the  subtraction  of  a  weighted  sum  of  yn_i  and 
yn_2  from  3C_  i,  so  that  P  would  become 

fnm   xn_x  -  ayn.i  +    byn_2  m   Xn.r[(1.  X)yn_1  +Ayn_,]   (32) 

a  ♦  b 

where     |A|  £  1,      The   operator   for   the    return   link    is    then  of   the 

form   (1-  A)    ♦     Az.  .0   for   the    new   system. 

The    relation    (*tween  x"   and  y  may    be   writ  as 

y    =    (1";;!1)    ii=22  k   I    sgn      [x    -    11- A*  Xx)y]  (33) 

s  £   i-  J 
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(1-A)*AZ 


Mgh 


71  r 

Z  H 


l-i 


1  -  zn 


>  8.   Proposed  dolta  i.-.odulat ion  for  ndnission  of  rnrrps 
using  delay  in  feedbaok  rotum. 
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A  recursion  relation  may  be  obtained  from  this  equation 
in  the  same  manner  that  equation  (14)  was  obtained  from 
equation  (7)  and  which  is 

Vn  ■  yn-l  *  m  {sSns[*n-l  "  (1"  ^n-l  "  A  yn-2] 

"  sSng[xn-m-l  "  (1-A>yn-m-l  "^n-m-2]|   <34> 
By  setting  A=  0,  one  may  calculate  the  unity  feedback  case. 

Various  values  of  A  between  0  and  1  may  then  be  tried  to  find 
conditions  for  optimum  performance  corresponding  to  best 
agreement  of  the  output  with  the  input. 

Computer  Results  for  Various  Combinations  of  Parameters 
with  (1-  \)+ X  z  Feedback.   For  the  particular  case  considered 
the  slope  of  x  was  set  equal  to  1/512T.   The  values  of  A  used 
were  0,  0.5,  and  1,  with  dead  spaces  in  the  signura  function  of 
0,  l/512 ,  and  l/l28  for  each  value  of  A  .   This  combination  thus 
included  the  unity  feedback  case  with  and  without  sipnum  dead 
spaces   and  covered  the  extremes  and  the  midpoint  for  A  in  the 
(1-  A)+ Az  compensation. 

See  rigs.  9-17,  Plates  I  through  III,  for  plots  of  the 
output  for  A=  0,  0.5,  and  1,  wirh6=  0.   The  output  fluctuated 
about  the  input  in  a  regular  "attern  with  abrupt  pattern  ch.. 
every  257th  sample.   Therefore  only  segments  of  each  consecutiv 
set  of  257  samples  were  plotted.   Analysis  of  the  outout  she 
that  in  .  .  ■u.-ral  the  average  error  remained  constant  dlirln     6h 
257  sample  section  and  iiu  om  one  1      >n  to 

vi.ition  of  the  error,  a  measure  of  the  noise,  also 
re;  tant  during  each       on  and  increased  substantial- 

ly from  one  ■<     >n  to  the  next  . 


/ION  Or    PLATE  I 

Delta  modulation  with  unity  feedback.   Output  for 

at ,  x  (t)  =  t/512  . 
Pig,      Output  between  samples  500  to  560. 
Fig.  10.   Output  between  samples  240  to  300. 
Fig.  11.   Output  between  samples  0  to  80. 
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PLATS  I 


Fig.  9. 


Fig.  10. 


Sample  numbor 
Fig.  11 


..  N  OF  PL     IX 

:a  modulation  with  (l-A)*Az  feedback,  where  ^  =  0.5. 
Output  for  input  x(t)  =  t/512T. 
Fig.  12.   Output  between  samples  500  to  560. 
Fig.  13.   Output  between  samples  240  to  300. 
Fig.  14.   Output  between  samples  0  to  80. 
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PLAT 2    II 


Fig.    12 


Fig.    13. 


Sar.  lb  or 

•  Lg«    111. 


Delta  modulation  with  (1-  A)+ Az  feedback  whore  \=    1.0. 

Output  for  input  x(t)  =  t/512T. 

Pig,  15.   Output  between  samples  500  to  500. 

Fig.  16.   Output  between  samples  240  and  520. 

Pig,  17.   Output  between  samples  0  to  80. 
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PLATS    III 


Fig.    15 


Fig.    16. 


Sample  nic.be r 
Fig.    17. 
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rage  error  was  calculated  by  a\     ing  the  difference 
between  x  and  y  at  each  of  the  sampling  ir.     Is  over  the 
period  of  repetition  of  the  fluctuations.   That  is 

(€>  -  £  E  cxi  -  yi>  ■  i  r  Si  (35) 

1=1         k  1=1 


The  standard  deviation  of  the  error  is 


<7-V*  £«/-£) 


~~1 

2 


(36) 


See  Table  3  for  a  tabulation  of  nver^.c  errors,  maximum  errors 
and  standard  deviations  of  the  error  for  A=  0,  0.5,  and  1, 
which  are  showTi  for  three  regular  sections  of  257  samples. 

Table  3  shows  that  the  effect  of  delay  in  the  feedback 
return  is  to  increase  the  noise,  the  maximum  error,  and  the 
average  error.   Reference  to  the  curves  of  output  vs.  input  with 
A^O  shows  that,  as  the  number  periods  increased,  the  output 
continued  to  attempt  to  reach  the  input  but  could  do  so  for 
only  one  or  two  sampling  intervals.   It  then  fell  behind  even 
further  than  for  the  unity  feedback  case  and  remained  behind 
for  a  considerable  time.   This  indicates  why  the  average  value 
of  the  output  decreased  for  larger  A . 

The  results  of  the  calculations  with  b=  tj-j-J  anc*  Y7E   snowed 
that  the  effect  of  8  was  to  smooth  out  fluctuations  at  the 
expense  of  increased  delay.   With  larger  values  of  A,  the 
fluctuations  increased  considerably  but  remained  loss  than 
those  for  6=0.   However  the  decreased  noise  did  not  seem  to 
be  worth  the  added  delay  except  for  the  case  with  6=  l/512. 
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Table  3.   Maximum  error,  average  error,  and  standard 
deviation  of  error  for  input,  x(t)  =  r/512, 
(1-A)+Az  in  feedback  link. 


A 


n 


Max.  -,rror 


Av£.  -3rror 


Std.  Dcv.  Error 


1-257    5.8594xl0~3 

0.0    258-514    9.7656xl0~3 

515-771   15.6250xl0"3 


0.9766x10 


-3 


0.9766X10--5 
10.2539xl0-3 


3.5210xl0"3 
6.5510xl0"3 
6.0001xl0"3 


1-257    7.8125xl0~3      0.9766xl0"3 

0.5    258-514   19.5313xl0"3      2.9297xl0~3 

515-771   23.4375xl0"3     12.6553xl0-3 


4.4752xl0"3 
9.8l43xl0~3 
9.8143X10"3 


1-257   13.6719xl0"3      0.9766xl0"3 

1.0    258-514   25.3906xl0"3      4.3S28X10"3 

515-771   37.1094xl0~3     21.8099xl0"3 


6.9283xlo~3 
14.77faxl0"3 
13.362xl0"3 


Table  4.   Maximum  error,  aver.'.      :-or,  and  standard 
deviation  of  error  for  input  x(t)  =  t/512. 

*ad  spaco   =  l/512  i;     num  function,  and 
unity  foodback. 


n 


. 


Ay.  .    -or 


ov.  ;rror 


1-257 
258-514 
515-771 


>xl0' 


,"3 
.  I393xl0"3 

7.8l25xl0"3 


.    I        '<Xl0 
2.^2^8x10 


.0x10 


-3 


-3 


-3 


2.18xl<  " 

.   .  |     x.O"3 

...      ■  .      " 
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ri   .   i-20f  Plate  IV,  show  the  output  for  S=  l/512,  which 
ovided  the  smoothest  output  with  tta     .st  standard  deviation 
and  least  maximum  error  for  any  case.   If  one  is  willing  to 
tolerate  an  average  error,  one  will  have  a  minimum  amount  of 
noise  throughout  the  thre»    |  iods  for  A=  0,  6=  1/512.   Calcula- 

;\s  show  that  the  standard  deviation  of  the  error  remained 
con      for  the  three  periods.   Sec  Table  4. 

^s  seems  to  be  as  nearly  ideal  as  can  be  obtained  using 
this  particular  arrangement.   The  conclusion  drawn  from  this 
set  of  calculations  is  that  the  basic  delta  modulation  system 
will  almost  follow  properly  chosen  ramps.   The  amount  of  delay 
between  input  and  output  is  small,  but  tends  to  increase  with 

ch  degenerate  triangular  wave  perio  .   The  amount  of  noise 
can  be  minimized  by  the  use  of  a  small  dead  space  in  the  signum 
function.   The  use  of  a  delay  in  the  feedback  link  produces 
no  beneficial  effects,  but  does  increase  delay  and  noise. 

Recursion  Formula  and  Results  for  l/(l+Ts)  in  Feedback  Link. 
The  next  compensation  scheme  employed  an  approximation 
to  the  identity  operator,  l/(l+Ts),  in  place  of  the  unity  in 
the  return  link.   This  procedure  is  analogous  to  the  compensation 
employed  in  continuous  servos  to  admit  ramp  inputs.   The 
resulting  system  is  shown  in  Pig*  21. 
\e.   equation  of  this  system  is 


_   q-z)(i-zm)    r 
y  =  ; k  Z  sgn 

mTs^  I 


x  - 


1  ♦  Ts_ 


(37) 


EXPLANATION  OF  PLATC  IV 

Delta  modulation  with  unity  feedback  and  signum  dead 
space  of  magnitude  1/512.   Output  for  input  x(t)  =  t/512T. 

Fig.  18.  Output  between  samples  500  to  580. 
Fig.  19.  Output  between  samples  240  to  320. 
Fig.  20.   Output  between  samples  0  to  80. 
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?LA 


Fig.    18. 


Fig.    19 


Sacplo  nur.ber 
lg.    20. 
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The   z -transform  of   both  sides   of  the   equation  is   taken  giving 


_       k   l-zra 
zy  =  m  TTT"     sSn 


zZx    -    zZ 


1+  Ts/ 


(38) 


Use    the   trapezoidal   approximation   (Kalijak    (5))    to   obtain 

Z(tg)    =  T(Zf)(Zg)    -   |T   [f(0)Zg    ♦    g(0)Zf]  (39) 


Consider 


Z(fg)    =   Z 


J^Ts  <*>] 


whe 


rein  f   =  y,    g  =         Ta    ,    f(0)    =  yQ   =   0,    and 


g(0)    =   L 
Also   needed   is 


-1 


1+  Ts 


oo 


=  I  *-t/<* 


t=0 


t  =  0 


1 
T 


(40) 


[l+TsJ       *      r=1 


(41) 


Hence 


1      y 


I*  Ts, 


_        oo 


=   T 


^X       (e-T/^z)r   Zy    -   i-  Zv 


L      r=l 


(42) 


Finally 

k    (l-zm^ 


-  oo 


zZx    -   ^.  z 


L  r=0 


"T/^r   .   1 


(e"1/^) 


2y 


(43) 


And 


yn  =  yn*l    *  H  sSn  [*n-i    ♦     5**11-1    "   |    Z     (c~T/r>ryn-r-l] 

*-  r=o  J 

oo 


m 


'n  [*n-m-l  •  ^  yn-m-l  "  F  Z  ^'T/T)ryn-r-m-l] 

*-  *■   r=o  -J 

(44) 
In  the  calculation  of  outputs,  it  was  not  known  what  vnlu* 
C  would  provide  an  optirvim  response,  so  the  following  range 
of  values  was  used:  X  ■  L92,  128,  64,  32,  16,  0,  4,  and  2.   Th* 
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LlM  of  6  was  0,  an        lope  of  the  input  was  l/512'f,  with 
ad  ra  unchanged. 

results  were  not  without  compromise.   They  showed  that 
output       can  be  eliminated  at  the       a  of  Lncn 

and  the  appear  Lay  during  the  later  sampling 

interval  .   The  output  pattern  was  similar  to  previous  cases 
in  that  it  formed  sections  of  257  samples  in  which  the  fluctua- 
tion remained  about  the  same,  although  not  as  regular  as  before. 
The  pattern  changed  from  one  section  to  the  next,  usually  with 
sharply  increasing  nois*  . 

The  optinim  value  of  T      ired  to  be  in  the  neighborhood 
of  A.   See  Figs.  23-25,  Plate  V.   The  fluctuation  became 

lent  for  larger  values  of  T,  and  delay  began  to  apj  ear  for 
T  loss  than  4.   The  minimum  fluctuation  for  any  case  calculated 
was  for  T=  2?  but  was  still  greater  than  the  corresponding 
viation  for  the  cases  with  (1-A)  *  Xz  compensation.   For 
example,  withT=  2,  n<257;  0^=  9.00  x  10" 5  while  the  corresponding 
deviation  for  the  worst  case  of  the  (1-X)  +  Xz  compensation, 

•  •re  A=  lf  n<257,  was  6.93  x  lO--*.   However  the  average  error 
for  both  cases  was  equal  to   0.9766  x  10"  .   Further  decreases 
in  T  would  likely  result  in  performance  approaching  that  of  the 
uncompensated  case. 

A  tabulation  of  tho  maximum  error,  average  error,  and 
standard  deviation  for  the  output   <  etiions  and  f or  'r=  A  is  given 
in  Table  5.   It  can  be  seen  from  the  tabulation  that  the  output 
tended  to  lead  the  input  for  the  first  257  samples,  and  then 
began  to  fall  behind,  although  not  significantly  until  n  became 
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Fig.  21.   Delta  modulation  with  y+rs~  operator  in  return  link. 


Table  5«     Klmum  error,  average  error,  and  standard  deviation 
of  error  for  input  x(t)  =  t/5l2T,  v/ith  l/l*rs  in 
feodback  return,  t=  ];. 


n 

;  !.'ax.    Error 

AV5.    Error 

':  Std.    Dov.    Error 

1   -   257 
296  -  3l\h 
566  -  62U 

21.U.   x  10"3 
Ul.O  x  10"3 
37.1   x  10"3 

-6.81;  x  10-3 
0.72  x  10"3 
8.86  x  10"3 

9.21  x  10-3 
17.58  x  10"^ 
15.5.0  *  10"-> 

1     +T3 


— 


1    -     X 

roTa 


■ 


Pig.    22.      Deltn  modulation  with  l*ra   Operator   foodlr  nun 

funotion. 


EXPLANATION  OF  PLAT£  V 

Delta  modulation  with  l/(l+Ts)  feedback  where  t  =  4 . 

Output  for  input  x(t)  =  t/512T. 

Fig.  23.   Output  between  samples  500  to  500. 

Fig.  24.   Output  bet    .  samples  240  to  320. 

Fig.    25,      Output    between  samples    0   to   80. 
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Fig.    2I4.. 
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gn  514.   Coir.parison  with  previous  tabulations  shows 

that  the  del.  lor  all  prcviou  :.  tor 

T"  =  4  except  tor  the  uncompensated  delta  modulation  with  a 
small  dead  space  in  the  si^num  function  and  n  >  514.   There 
the  avc i     error  was  only  4.88  x  10'5.   For  the  first  two 
output  sections  the  error  for  the  l/(l*4s)  case  was  better.   With 

d  to  considerations  of  £   and  07  ,  the  uncompensated  case 
was  better.   Its  maximum  error   as  l/3  that  of  the  l/(l«-4s) 
case,  and  its  noise  was  17  db  under  that  of  the  l/(l+4s)  case 
for  n  >  514. 

If  the  preat  increase  in  noise  can  be  allowed,  an  ov< rail 
compromise  on  the  delay  can  be  reached  by  proper  choice  of  T" , 
perhaps  yielding  a  zero  overall  average  delay.   The  amount  of 
noise  that  can  be  tolerated  and  the  final  compromise  made  will 
depend  on  the  application  at  hand. 

Derivative  Compensation  of  Jrror.   Another  method  of 

■         ■      i        i   -         i 

compensation,  insertion  of  a  function  of  the  form  (1+ts) 
between  the  comparator  and  the  signum  function,  was  attempted 
next.   Set'   .  .22.   The  equation  of  this  system  is 

(i-zn)  /i-z\     r      _  _n 

y    =   m       Ts  [—)  Z   s^    [<l*T«)(x-y)J  <45) 

The    solution  of   this   equation  yields    the    following 
recursion   formu. 

*H  =   *U-1    *  ■  s*n  [<T*  *    1)(xn-l    -  yn-l> 

oo 

-  ^T     rT2    <-1>r^^-r->'n-r)]    -   B   s^n  [^  +1  >(xn-m-l 

"   *n-m-l>    "  f     ZL    ("1>r   ^-m-r^n-m-P] 

r=2  -> 
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Again,  proper  values  for  t-  were  unknown,  and  values  of  16, 
8,  4,  2,  1,  -,  i,  — ,  and  ~  were  used.   Slopes  of  1/512T  and 
1/256T  for  x  were  used,  and  S  was  set  equal  to  zero.   Other 
parameters  remained  unchanged . 

The  effect  of  this  modification  was  to  introduce  consider- 
able delay,  increasing  with  T,  with  a  certain  amount  of 
fluctuacion.   For  t">l/4 ,  in  every  case,  the  output  '..-as  very 
small,  hovering  slightly  above  zero,  since  the  error  was  almost 
equal  to  the  input.   It  could  be  said  that  th^  output  consisted 
of  the  noise  generated  plus  a  small  constant.   For  T'a  l/8,  l/l6, 
the-  output  resembled  a  ramp  with  the  error  remaining  more 
nearly  constant.   For  T"=  l/8  the  average  delay  was  approximately 
175T  while  for  T  =  l/l6  it  was  approximately  75T.   This  much 
delay  is  intolerable  and  the  fluctuation  was  still  greater  than 
is  desired.   It  was  concluded  that  this  particular  modification 
was  unsuccessful. 

Conclusions  on  Admission  of  Ranps.   The  results  of  the 



investigation  of  delta  modulation  response  to  ramp  inputs  led 
to  the  following  conclusion.   A  trade-off  exists  between  error 

.  noise,  and  an  ideal  outout  of  zero  error  and  constant 
minimized  fluctuation  is  unattainable.   If  one  desires  minimal 
noise  and  is  willing  to  tolerate  small  errors,  a  small  dead 

ce  in  the  signum  function  of  a  deli  on  system  with 

unity  feedback  will  provide  such  performance.   On  the  Othi     ad, 
if  zero  ror  is  desired  with  li  m  for  flu 

tion,  then  such  pcrfon        n  i*-  ob:     d  through  the  u.se 
l/(l*  Tj;)  fcedb      1th  T   prop.        >>s.n. 


43 


Review  of  the  Goals  and  Capabilities  of 
...  •  ic  and  I>cterministic  An     s 

deterministic  analysis  of  the  delta  modulation  system 
has  been  presented.   It  would  be  of  interest  to  compare  the 
accomplishments  and  capabilities  of  this  approach  with  those  of 
the  probabilistic  analysis  of  Zetterberg.   Since  do  Jager's 
work  was  primarily  descriptive  an        imcntal,  there  is  no 
n-^ed  to  discuss  it  further. 

The  goal  of  the  probabilistic  analysis  was  to  determine  the 
operties  of  delta  modulation  as  an  information  encoding- 
decoding  system  from  the  viewpoint  of  information  th.  ory  and 
to  compare  these  properties  to  corresponding  ones  for  pulse  code 
modulation.   Considerable  emphasis  was  placed  on  this  comparison 
rather  than  evaluating  DM  as  a  new  system.   Among  the  properties 
considered  were  channel  capacity,  signal  power,  optimum  ensemble 
and  its  spoctral  properties,  and  noise  power  generated  within 
the  system. 

The  probabilistic  approach  has  th.  capability  of  giving  an 
iication  of  how  well  a  communication  system  will  oerform  with 
delta  modulation  as  an  integral  part.   This  is  carried  out  by 
computing  the  amount  of  information  that  can  be  handled,  given 
requirements  regarding  quality  of  demodulated  signals.   The 

itations  on  types  of  inputs  that  can  be  reproduced  with  desired 
accuracy  and  the  optimal  inputs,  together  with  upper  bounds  on 

it  frequencies,  may  be  obtained.   An  idea  of  channel  bandwidth 
requirement  may  also  be  computed.   By  this  means,  one  can  discover 
the  types  of  applications  for  which  DM  is  best  suited  along  with 
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ideas  on  its  incorporation  into  larger  systems.   The 
probabilistic  analysis  deals  with  the  overall  performance 
of  DM  and  is  incapable  of  predicting  specific  events  at 
particular  times.   It  is  unable  to  determine  least  upper  bounds 
on  information  capacities  leaving  the  user  uncertain  as  to  how 
far  away  actual  capacities  are. 

This  uncertainty  points  out  the  need  for  other  methods 
of  studying  the  system,  and  the  deterministic  analysis  of  this 
paper  fills  a  part  of  this  need.   The  goal  was  to  obtain  the 
magnitude  of  the  demodulated  output  as  a  function  of  time, 
given  the  input  as  a  function  of  time,  and  to  provide  compensa- 
tion schemes.   Deterministic  analysis  is  capable  of  dealing  with 
absolute  quantities  at  suecific  times,  and  furthermore  a  measure 
of  the  overall  performance  as  well  as  performance  in  a  larger 
system  may  be  derived.   Also  problems  encountered  in  the 
physical  design  of  a  delta  modulation  system  require  deterministic 
methods  for  adequate  solution. 

The  exact  form  of  the  output  of  the  demodulator  can  be 
calculated.   From  thlfl  one  can  eacerl  ain  the  faithfulness  of 
reproduction  along  with  signal  to  noise  ratios.   Since  the  exact 
form  of  the  modulator  output  is  also  uv.iilabl  ,      :.el 
requirements  and  bandwidths  may  te  obtained.   The  di-it 
computer  r..     it  possible  to       out  the  arithmetic  opera- 
tions o         term;.     .  c  method  very  •      •  ,  ;,  ny 
time  element  drawback.   This  method  also  p     <-a   an    I  llent 
as  for  evaluating  the  effect!  of  new  modifications  on  I  I 
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i-tn,  as  is  ill      red  by  servo  comp.     ion  t<       >:s 
cribed  earli*-r.   In  most  cases   a  probabilistic  analysis 
woul<  o  information  of  this  sort.   Although  t: 

fu.     as  of  the  two  approaches  overlap  slightly,  the  fact 
remains  that  both  are  necessary  for  a  complete  pictu:  . 

[ON  OF  DSLTA  MODULATION 

It  is  possible  to  generalize  a  delta  modulation  system  by 
replacing  the  three-level  quantizer  with  a  five-level   quantizer. 
For  basic  development  of  the  new  system   reconsider  DM  as  the 
original  combination  of  pulse  amplitude  and  pulse  duration 
modulators  and  demodulators.   Let  the  input  signal  be  fed  into 
two  separate  PDAM  modulators.   Both  of  these  modulators  consist 
of  the  usual  comparator  with  its  triangular  wave,  followed  by 
a  signum  function  and  a  sampling  operator.   For  section  1 
rrodulator  the  signum  function  has  no  dead  space,  but  for  section 
2  it  has  a  dead  space  of  magnitude  D.   The  triangular  waves 

■  both  of  the  same  period  bat  have  dlf     >t  amplitudes.   Now, 
let  the  two-pulse  duration-amplitude  modulated  signals  be  added 
to  yield  a  five-level  sequence  of  pulses.   This  combination  of 
PDAM  modulators  operating  in  tandem  makes  up  the  new  modulator. 

Since  the  PDAM  demodulators  for  each  section  are  identical, 
the  use  of  a  single  demodulator  for  the  added  outputs  is 
equivalent  to  demodulating  each  section  separately  and  adding 
the  resultants,  and  this  justifies  addition  of  the  two  signals. 
A  diagram  of  the  new  arrangement  is  shown  in  Fig.  26. 
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An  analysis  of  each  section  rr.ust  be  made  in  order  to 
obtain  a  combined  input-output  characteristic  which  is  linear. 
For  simplicity  it  will  be  assumed  that  the  PAM  portion 
reproduces  perfectly  so  that  only  the  PDH  part  need  be  analyzed. 
Since  the  modulator  in  section  1  has  no  dead  space  in  the 
signum  function,  when  combined  with  the  demodulator  it  gives 
an  overall  input-output  characteristic  which  is  linear,  passes 
through  the  origin,  and  saturates  at  output  level  mkT,  where  mT 
is  the  period  of  the  triangular  wave  and  k  is  the  amplification 
factor.   A  plot  of  this  characteristic  is  shown  in  Fig.  27. 
Saturation  occurs  for  an  input  magnitude  equal  to  H^,  where  H^ 
is  the  amplitude  of  the  triangular  wave  of  section  1.   Hence  the 
slope  of  the  output  versus  input  curve  is  equal  to  mlcT/H^.   It 
must  be  pointed  out  that  these  response  characteristics  apply 
only  to  stop  inpuus. 

It  is  desirable  that  the  response  characteristic  of 
section  2  be  such  that  it  will  complement  the  response  of 
section  1  to  produce  a  linear  resultant  when  the  two  outputs 
are  added.   Figure  28  shows  the  desired  output  versus  input 
curve  for  section  2.   It  has  a  dead  space  for  which  the  output 
is  zero  lor  inputs  of  magnitude  less  than  some  value , A .   For 
inputs  greater  than  A  the  output  rises  linearly  until  an 
output  saturation  level  is  reached.   An  at  ilar  to  the 

one  for  section  1  shows  that  such  a  char     ist  ic  mav  be 
obtained  only  when  the  magnitude  of  the  dead  space ,   ,  of  the 
signum  function  is  greater  than  the  amplitude  o  alar 

vc,  H2.   The  value  of  A  ,  which  I  .  x,  :or 
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Pig.  26.   Dual  section  pulse  duration-aj-plituic  modulation. 
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Pig.  27.   Response  of  PL       Fig.  28. 
•.vie.. out  do  ad  spaco. 
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Pig.  29.   Response  of  dual  section  PDM. 
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which  output  begins,  equals  the  difference  between  the  magnitude 
of  signum  dead  space  and  the  triangular  wave  amplitude,  or 
A  -   D-H2.   Saturation  occurs  for  an  input  equal  to  (D+H2).   The 
output  saturation  level  is  again  mkT  so  that  the  slope  of  the 
out put -vers us -input  line  is  mkT/2H2.   In  order  to  obtain  a 
linear  overall  characteristic,  D  and  H2  must  be  chosen  so  that 
section  2  output  begins  where  section  1  output  saturates,  which 
implies  that  (D-H2)  must  equal  H,.   Furthermore  the  slopes  of 
the  two  characteristics  must  be  equal,  or  mkT/211^  =  mkT/2H2# 
Therefore  K2  =  Hly/2 ,  and  D  =  5U1/2.      The  overall  transfer 
characteristic  is  shown  in  Fig.  29.   This  combined  system  will 
now  respond  to  inputs  whose  magnitudes  range  from  zero  to  (D+H2) 
2H^,  with  a  maximum  output  equal  to  2mkT.   The  effect  has  been 
to  double  the  range  of  response  of  pulse  amplitude  duration 
modulation,  while  retaining  the  same  accuracy  of  response,  since 
the  pulse  duration  modulator  period,  mT,  and  pulse  amplitude 
modulator  period  T,  have  remained  unchanged.   If  th        ampli- 
tude range  were  maintained  the  accuracy  would  double. 

Now  let  the  triangular  waves  be  removed  from  both  sections 
of  the  modulator,  and  let  the  output  of  the  demodulator  be 
substituted.   Since  the  amplitude  of  the  triangular  wave  for 
section  2  is  to  be  l/3  the  magnitude  ot  the  dead  space,  n 
limiting  amplifier  may  be  placed  in  the  return  to  thi     I  tion 
which  will  limit  the  amount  Ted  back  to  a  : 

of  D/3.   This  ik-w  0O1  own  in  P   .  50.   An  ev. 

furth. ■:•  •...  .  . •:  ing  of  the  system  would  be  t 
limiting  amplifier.   -  Ll  than  woul       Lt  sir        lion  I 


49 


~&- 


+  ii 


-  ZtfSH 


l 

-    z 
a 

k 

1  -  zn 


-♦ — *- 


.  30.   I>ual  section  PDA?.'  weakened  to  a  feedback  dovico. 
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the  new  system  to  the  old  delta  modulation  with  the  signum 
function  replaced  by  a  new  five-level  quantizer.   The  quantizer 
is  the  result  of  adding  the  outputs  of  the  ordinary  signum 
function  and  the  signum  function  with  a  dead  space.   See  rig.  31. 

This  weakened  FDAM  with  its  multilevel  pulses  is  equivalent 
to  a  digital  servo,   however,  accurate  demodulation  after 
passage  through  an  imperfect  transmission  channel  requires 
regeneration  of  pulses  to  unit  height.   This  can  be  achieved 
by  representing  a  pulse  of  the  higher  level  by  two  pulses  of  a 
single  level,  with  the  second  pulse  occurring  halfway  between 
the  original  pulse  and  its  following  pulse.   Hence,  if  the  pulse 
to  be  encoded  occurs  at  time  t  =  nT,  where  n=  1,  2,  3,  ..., 
the  first  pulse  of  its  coded  representation  occurs  coincidently 
at  t  =  nT,  and  the  second  at  t  =  nT  ♦  T/2.   The  same  situation 
will  apply  to  negative  pulses. 

n  encoder  for  this  .urpose  is  shown  ir.     .32.     e 
five-level  sequence  of  pulses  is  fed  into  the  device  and  passes 

nura  function  with  a  dead  space  such  that  it  will 
produce  an  output  for  only  the  higher  l«-vel  posit i  ;Lve 

Iscs.  The  output  of  the  signum  function  is  then  di 
time  t  =  T/2  and  added  to  the  original  input  sequence.   This 
sum  is  passed  through  a  limiting  amplifier  (•       nt  to 
another  signum  function)  to  obtain  the  l  in.il  encoded  output 
sequence  of  pulses  with  m  uniform  m.         (and  null   ulses 
also) . 

The  decoder  is  even  simpler.  it  it  will 

a  sequence  ot  null  or       rm  mngnitu  ora 
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.  33*   Decoder  for  encoded  five-level  delta  r.odulation. 
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suitable  channel  amplifiers.   This  sequence  is  delayed  by 
time  t  =  T/2  and  added  to  the  undelayed  sequence.   This  sum 
is  then  amplitude  sampled  at  times  t  =  nT  ♦  T/2,  (n=l ,2 , 3 , . . . ) , 
so  that  the  original  sequence  of  five-level  pulses  is 
reconstructed,  although  delayed  by  an  amount  T/2.   See  Fig.  33. 

It  has  been  shown  that  a  five-level  delta  modulation  is 
feasible,  and  that  the  sequence  of  pulses  comprising  the 
modulated  signal  can  be  constrained  to  have  a  single  magnitude, 
which  retains  a  desirable  feature  of  DM.     •  accuracy  of 
reproduction  has  been  doubled  for  a  given  ran^e  of  inputs.   The 
price  for  this  improvement  is  doubled  bandwidth  required  for 
a  given  fundamental  sampling  frequency,  or  a  hplving  of  the 
fundamental  sampling  rate  for  a  given  bandwidth  restriction  of 

•  channel.   This  is  caused  by  insertion  of  the  extra  pulses 
in  the  coded  output.   An  increased  delay  in  the  output  has  also 
been  incurred.   i-\irther  increases  in  the  number  of  signum 
levels  will  lead  to  a  pulse  code  modulation  for  encoding  the 
delta  modulated  signal. 

Analysis  of  the  five-level  delta  modulation  to  determine 
how  far  the  PDAM  may  be  d>         id  the  best  c'r.oicc  of 

ametc:     s  not  been  carried  out  Lm*«   However,  an 

idea  of  the  performance  of  a  ■ultilftVftl  system  has  been  prov 
by  David  (3). 
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summary 

...tory  of  previous  work  in  the  field  of  delta  modulation 
has  been  presented  along  with  an  interpretation  of  results. 
Past  writers  considered  DM  to  be  similar  to  pulse  code 
modulation  and  used  the  latter  as  a  yardstick  to  measure  the 
performance  of  DM.   The  analyses  were  probabilistic  in  nature. 

upper  bound  on  channel  capacities  was  derived  as  well  as 
noise  and  spectral  calculations. 

Barber  interpreted  delta  modulation  as  a  combination  of 
pulse  duration  and  pulse  amplitude  modulation.   This  provided 
the  basis  for  a  deterministic  analysis  of  DM  lifting  Z-transform 

nods  and  digital  computer  techniques.   The  results  of  these 
Calculation*  showed  the  system  to  have  the  characteristics  of 
a  servo  for  certain  ranges  of  step  inputs.   Signal-to-noise 
ratios  were  calculated  from  the  computor  results.   Also  the 

ct  form  of  the  modulator  output,  a  sequence  of  positiv  , 

alive,  and  null  pulses,  was  comnutc  . 

An  investigation  of  th  ncc  of  DM  for  ramp  inputs 

was  made,  and  the  effects  of  compensation  techniques  were 
studied.   It  was  found  that  for  DM  with  unity  feedback,  the 

st  satisfactory  performance  could  be  obtained  through  the  use 
OL  a  signum  function  with  a  small  dead  space.   This  minimized 
output  fluctuations  although  a  small   increasinf  error  remained. 
3y  placing  a  function  of  the  form  l/(l*ts)  in  the  return  li. 

verage  error  could  be  eliminated  at  the  expense  of  rro at ly 
increased  noise.   Further  investigation  slowed  that  an  idealized 
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resoonse  of  minimal  noise  and  zero  error  is  unattainable  so 
that  a  compromise  will  have  to  be  made  between  the  two. 

A  comparison  of  the  goals  and  capabilities  of  the 
probabilistic  and  deterministic  analyses  showed  that  both  were 
necessary  for  a  more  nearly  complete  picture.   The  deterministic 
approach  was  able  to  supoly  specific,  concrete  results  and 
was  best  suited  for  use  in  analysis  and  design.   The 
probabilistic  approach  was  able  only  to  set  upper  bounds  on 
quantities  and  was  better  suited  for  evaluating  DM  as  an 
information  handling  system. 

Finally,  an  extension  of  delta  modulation  using  fivv- -level 
quantization  was  presented.   The  new  system  was  developed  from 
the  rudimentary  combination  of  pulse  duration  and  pulse  amplit; 
modulation.   It  should  reproduce  the  input  more  accurately 
over  the  same  range  of  input  values.   An  encoder  :       d  to 
change  the  five-level  pulses  into  positive  and  I     Ivc  puli- 
of  a  single  magnitude,  effecting  a      -alization  of  delta 
modulation. 
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Previous  investigators  have  made  probabilistic  comparisons 
of  delta  modulation  with  pulse  code  modulation.   However  Barber 
has  shown  that  delta  modulation  is  actually  a  hybrid  of  pulse 
duration  and  pulse  amplitude  modulation.   Let  a  PDM  modulator 
and  its  demodulator  be  separated  by  a  PAK  modulator  and 
demodulator.   Then  remove  the  triangular  wave  from  the  pulse 
duration  modulator  and  replace  it  by  the  demodulated  output. 
This  weakened  device  is  a  delta  modulation  system. 

An  idealized  mathematical  model  may  Le  derived  and  leads 
to  a  deterministic  analysis.   The  weakened  PDM  modulator  is 
represented  by  a  comparator  followed  by  a  signum  function.   Its 
demodulator  is  a  holding  filter  whose  transfer  function  is 
(l-zm)/mTs.   The  PAM  modulator  is  the  sampling  operator,  Z.   Its 
demodulator  is  another  holding  filter  with  the  transfer  function, 
(l-z)/s,  which  functions  as  a  linear-phase  low-pass  filter.   A 
Z-transform  analysis  leads  to  a  recursion  relation  which  enabl. 
one  to  compute  the  demodulator  output  at  each  sampling  instant , 
given  the  input  as  a  function  of  time.   Digital  computer 
techniques  applied  to  the  recursion  fromula  show  that  the  basic 
delta  modulation  system  can  be  unction  as  a  servo  for 

a  range  of  step  inputs.   PrOB  these  results  o     .y  compu 
signal-to-noise  ratios.   The  exact  form  of  the  modu1.      signal 
may  also  be  computed. 

An  e      I  on  of  the  iu-       tion  co;.  vsponse 

of  a  DM  system  to  ramp  inputs  and  provides  proper  mo        on 
for  fcl  .  ra;..:   .   It  is  found  t 


Ill  almost  follow  ra  .  cons.         .nount 

noise  and  a  small  error,  both  of  which  in<  iodically 

by  jv.        .     The  insertion  of  a  small  dead  space  into  the  signura 
function  will  minimize  the  noise  to  a  cc     ..t  value  althou 

8  error  incr*.      slightly.   A  normalized  delay  in  the  return 
link  to  eliminate  the  error  has  the  opposite  effect  of  increasing 
bo:     .or  and  noise.   An  Approximation  to  the  identity  operator, 
l/(l  ♦  Ts),  in  the  return  link  can  be  made  to  -rovide  zero 

rage  output  error  over  the  interval  of  interest  through  the 
proper  choice  of  T .   The  price  for  this  is  increased  noise. 

Further  study  shows  that  a  trade  off  exists  between  error 
and  noise  and  that  both  may  not  be  minimized  at  the  same  time. 
For  minimal  error,  use  l/( 1  ♦ Ts )  in  the  feedback  link,  and  for 
minimal  noise  use  unity  feedback  with  a  small  signum  dead  space. 

The  probabilistic  and  deterministic  analyses  are  both 
needed  for  a  complete  study  of  a  delta  modulation  system.   The 
former  provides  a  measure  of  the  capacity  of  the  system  to  handle 
ion.  The  deterministic  approach  is  necessary  for 
cific  analysis  and  the  design  of  a  system. 

Delta  modulation  may  be  generalized  to  a  form  having  five 
level  quantization.   This  should  result  in  improved  and  extended 

formancc.   Furthermore  it  is  possible  to  encode  the  five- 
lcvel  delta  modulated  signal  into  one  having  three-level 
quantization  by  representing  the  highest  level  pulses  by  two 
lower  level  nuls'-s.   This,  however,  results  in  doubled  channel 
bandwidth  requirements. 


